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Abstract The inverse of the relative increment rate of size 
is proposed as a measure of the maturity of a forest tree. 
The heartwood content within bole cross sections is mod¬ 
eled semiempirically using asymptotic fitting. It is shown 
that the present proposal for a definition of maturity statis¬ 
tically explains heartwood content far better than age or 
size. The heartwood content appears to be independent on 
growth rate and tree size. 
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Introduction 

In many tree species the permeability of heartwood to liq¬ 
uid and gas is less than that of sapwood . 1,2 In some tree 
species, heartwood has chemical constituents that inhibit 
fungal activity . 3 ' 5 These two properties make heartwood a 
material with natural decay resistance. Heartwood, with a 
significant extractives content, also has more stable dimen¬ 
sions when subjected to varying air humidity , 4,6,7 

Heartwood formation appears to take place only in trees 
exceeding a specified age; thereafter the proportion of 
heartwood increases with age . 68 ' 11 On the other hand, the 
amount of heartwood appears to be negatively correlated 
with the growth rate of any tree . 8 ' 611 ’’ 2 However, fast-grown 
trees have a lower heartwood content only when they are 
compared with more slowly grown trees of the same size . 13,14 
When the age is the same, large diameter trees appear to 
have a greater relative proportion of heartwood . 13-14 

The maturation of any plant is manifested in a variety of 
physiological changes. The sexual maturation of trees ap- 
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pears to be related to the accumulated number of cell divi¬ 
sions in the meristems, rather than age or size per se . 15-19 
Some observations suggest that pine cambium maturity, 
judged on the basis of the change in wood basic density, is 
better described in terms of chronological age than in terms 
of the. accumulated amount of production . 29 

It is well known that the decline in the relative increment 
rate of size due to age takes longer as growing conditions 
worsen. In other words, the relative growth at forest sites on 
poor soil or in a cool climate declines at an older age than 
the relative growth at sites where the soil is good and the 
climate is warm . 21 " 23 On the other hand, with favorable con¬ 
ditions for growth, trees tend to grow larger . 13,21-23 The rela¬ 
tive proportions of biomass compartments (heartwood, 
sapwood, foliage, branches, roots) change along with tree 
maturation in a similar pattern at good and poor sites, al¬ 
though at good sites such transitions take place at a younger 
age and larger size . 24 Thus, it appears that within a particu¬ 
lar tree species maturity cannot be described solely in terms 
of tree age, cambium age, or tree size. 

No tree grows to infinite size. Thus, at some stage the 
increment rate of size declines and ultimately approaches 
zero. Observations propose that the relative increment rate 
of size generally decreases with maturation . 25-27 Thus, the 
inverse of the relative increment rate of size generally in¬ 
creases monotonically. We propose that the inverse of the 
relative increment rate of size can be used as a measure of 
maturity. 

This proposed measure of maturity has the dimension of 
time. The ultimate outcome of maturation being death, the 
ultimate relative increment rate of size corresponds to zero 
per time unit, maturity thus approaching infinite time. On 
the other hand, because the size of a very young individual 
is small, the relative increment rate of size may be large and 
the proposed measure of maturity has a low value. The 
maturity of a newly established individual, however, does 
not correspond to zero. Also, this asymptote has a physi¬ 
ological correspondence: The establishment of an indi¬ 
vidual does not correspond to the origin of life. 

Heartwood formation is a phase transition, where tissue 
that includes living cells becomes a tissue that no longer 
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includes cells with any metabolism. Furthermore, consider¬ 
ing that the extreme outcome of maturation is death - a 
state where there is no metabolism - it is obvious that the 
extreme outcome of maturation refers to heartwood con¬ 
tent of unity. 

We intended to clarify experimentally to what degree 
heartwood content can be explained in terms of the pro¬ 
posed measure of tree maturity (i.e., the inverse relative 
increment rate of size). Heartwood content is modeled 
semiempirically using asymptotic fitting. The coefficient of 
determination by the achieved model is compared with the 
corresponding coefficients of determination by models 
based on age or size. We also intended to clarify whether 
the heartwood content at a specified maturity depends on 
the growth rate or the tree size. 


Materials and methods 

The experimental material consists of two independent 
datasets. Two datasets were collected to clarify the repeat¬ 
ability of the results. 

Dataset 1 was acquired during March 2000. One hundred 
Scots pine (Pinus silvestris L.) butt sawlogs entering a 
sawmill in eastern Finland were sampled. Sampling was 
designed to represent the volume-weighted diameter distri¬ 
bution and geographic yield of Kainuu Province. Sample 
disks were sawn 55 cm from the top of the sampled butt log, 
the logs being 3-6 m in length. Circular specimens (28mm 
diameter) were cut from the middle of the sapwood zone of 
the disk along a sampling line that corresponds to the great¬ 
est diameter of the disk. The average annual ring width 
within each specimen was measured, as was the distance 
from the tree pith, in terms of millimeters as well as in terms 
of the number of annual rings. 

Dataset 2 was acquired during August 2000. Forty logs 
were sampled from the same sawmill as was used to obtain 
Dataset 1. These logs also represented the volume-weighted 
diameter distribution and geographic yield of Kainuu Prov¬ 
ince. A sample disk was sawn from these logs at 90% of the 
log length, and specimens were produced as for Dataset 1. 

For the sample disks a coordinate system was established 
by taking the line of greatest diameter as the reference line 
and the center point of this line as the origin. The distance 
of the heartwood-sapwood boundary from the origin, as 
well as the distance of the wood-bark boundary, was mea¬ 
sured at angular intervals of n/6. The location of the geo¬ 
metric center of the heartwood cross-sectional area and the 
center of the log cross-sectional area were calculated; the 
location of the pith of the tree was also determined. 
The geometric area centers were taken where the sum of 
the squared differences of distances from the center to any 
of the boundaries was least; thus the geometric center did 
not exactly coincide with the pith. The cross-sectional area 
of the wood disks and the cross-sectional area of the heart- 
wood were determined on the basis of the average distance 
of the 12 boundaries from the center of the disk and from 
the center of the heartwood, respectively. 


Results 

Heartwood formation as a maturation phenomenon 

We have proposed the inverse of the relative increment rate 
of size as a measure of maturity. Such a quantity may 
achieve values between zero and infinity, having the dimen¬ 
sion of time. To ensure that such a quantity is readily mea¬ 
surable from any cross section of a bole, we implemented 
the maturity characterization in terms of the inverse of the 
relative basal area increment rate. Assuming that the cross 
section of any bole is essentially circular, this measure can 
be written as 


,« ■ trr 2 r 

Maturity = - z— = —— 

2 «r* 2*1 0) 

dr dr 

where r refers to the bole radius, and t refers to time. 

For the present experimental material, the maturity ac¬ 
cording to Eq. (1) can be determined from wood specimens 
obtained to represent the sapwood of the sample disk. This 
is because the annual ring width, which is needed for Eq. 
(1), was measured in the circular specimens cut from the 
middle of the sapwood section. The calculated maturity 
thus corresponds to the maturity that prevailed some years 
before felling, during formation of the xylem. which pres¬ 
ently constitutes the sapwood. The measured heartwood 
content, however, corresponds to the time of felling. Thus, 
there is a phase lag between the measurement of maturity 
and the heartwood content. The phase lag, however, is as¬ 
sumed here to be unimportant. 

The heartwood content of any butt log top cross section 
as a function of the maturity of the cross section is shown in 
Fig. 1. We find that despite the phase lag in the measure¬ 
ments the heartwood content is accurately determined by 
maturity. Any function describing heartwood content as a 
function of maturity must approach zero at a small finite 
value of maturity and approach unity at infinite maturity. 
The function 


Heartwood content = 0, if maturity < m 0 

/ 

Heartwood content = 1 — exp 

v 

if maturity 2 m 0 


maturity - m c 


( 2 ) 


where m 0 is the maturity where heartwood formation ini¬ 
tiates, and a and b are constants, satisfies these logical limit 
values. When fitted to Dataset 1, Eq. (2) provides a coeffi¬ 
cient of determination of 0.76 (Fig. 1). Dataset 2 displays a 
trend of the heartwood content increment as a function of 
maturity, similar to that of Dataset 1. 

Because maturity was defined above as the inverse of the 
relative increment rate of size, the heartwood content at a 
specified maturity does not need to be related to the width 
of the annual rings. Figure 2 shows that the normalized 
residuals of Eq. (2) do not correlate with the sapwood an- 
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Fig. 1. Heartwood proportion of a cross-sectional area in a wood disk 
taken from the vicinity of the upper end of a butt log as a function 
of maturity, evaluated using Eq. (1). Eq. (2) is fitted to Dataset 1 
(diamonds). Variation within Dataset 2 ( triangles ) can be reasonably 
described with the same model 



Fig. 2. Normalized residual variation of the heartwood proportion as a 
function of sapwood annual ring width after fitting Eq. (2) to Dataset 1 
(diamonds). At a specified maturity, the heartwood content does not 
correlate with the annual ring width. Variation in Dataset 2 ( triangles ) 
can be reasonably described with the same model as used for Dataset 
1 



Fig. 3. Heartwood proportion of a cross-sectional area in a wood disk 
taken from the vicinity of the upper end of a butt log as a function of 
the number of annual rings in the disk. An exponential function is fitted 
to Dataset 1 (diamonds). Variation within dataset 2 (triangles) can be 
reasonably described with the same model 


Dataset 1, provides a coefficient of determination of 0.46. 
Dataset 2 displays still more scatter. 

Figure 4 shows that once heartwood content has been 
predicted in terms of cambium age, an apparent effect of 
growth rate remains, the normalized residuals depending on 
annual ring width within the sapwood. The correlation with 
annual ring width is negative, which agrees with findings in 
the literature. 8,6 ' 11 ’ 12 We also find from Fig. 4 that at a speci¬ 
fied age the heartwood content correlates positively with 
size, which also agrees with the literature. 13 ' 14 

Heartwood content as a function of size 

Heartwood content as a function of log diameter is shown in 
Fig. 5. We find that size explains heartwood content still less 
accurately than age. An exponential function of the same 
type as Eq. (2), fitted to Dataset 1, provides a coefficient of 
determination of 0.31. Dataset 2 displays at least the same 
amount of scatter. Figure 6 shows that once the effect of size 
is removed, the heartwood content has a strong negative 
correlation with the annual ring width and a clear positive 
correlation with age. 


nual ring width, nor does the residual variation correlate 
with age or size. 

Heartwood content as a function of age 

Heartwood content as a function of the number of annual 
rings in a sample disk is shown in Fig. 3. We find that 
cambium age explains the heartwood content less accu¬ 
rately than does the maturity defined in Eq. (1). An expo¬ 
nential function of the same type as Eq. (2), fitted to 


Effect of growth rate on heartwood content 

Even if heartwood content at a specified maturity does not 
correlate with the width of annual rings (Fig. 2), it is un¬ 
known to what degree heartwood content possibly depends 
on the growth rate. This is because the absolute increment 
rate of size generally is not independent of maturity. 28 ^ 32 
Thus, the eventual independent effect of growth rate on a 
particular property of a tree or a tissue formed by the tree 
must be studied in terms of a dimensionless, maturity- 
independent measure of the growth rate. 27 
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Fig. 4. Normalized residual variation of heartwood proportion as a 
function of the sapwood annual ring width (top) and as a function of 
log diameter (bottom) 

A dimensionless, maturity-independent growth rate can 
be produced by normalizing an absolute increment rate of 
size by the increment rate of size that is statistically typical 
for that particular maturity. 20,27 In the present case, with 
maturity being defined as in Eq. (1), the growth rate is: 



Fig. 5. Heartwood proportion of a cross-sectional area in a wood disk 
taken from the vicinity of the upper end of a butt log as a function of 
the diameter of the disk. An exponential function is fitted to Dataset 1 
(diamonds). Variation within Dataset 2 ( triangles ) can be reasonably 
described with the same model 

the right-hand side of Eq. (3). Once maturity approaches 
zero, the bole radius also approaches zero. On the other 
hand, once maturity approaches infinity, the bole radius 
approaches a finite limit. Such a function is represented by 


I'm = c 

1 - exp 

( maturity 7 )' 'j 


_ 

v V d J 


where c, d, and e are constants, the constant c having the 
physical interpretation of the size limit. Eq. (4), fitted to 
Dataset 1, is shown in Fig. 7. Observations from Datasets 1 
and 2 are also plotted in Fig. 7. On average, the growth rate 
in Fig. 7 is somewhat higher in Dataset 2 than in Dataset 1. 

The effect of growth rate, as defined in Eq. (3), on the 
heartwood content is shown in Fig. 8. The heartwood con¬ 
tent is not dependent on the growth rate. 


Discussion 


Growth rate 



r 

r M 


( 3 ) 


where (dr/df) M is a typical (mean) value for the incremental 
rate of radial growth at maturity M, and is the typical 
(mean) bole radius at maturity M. The equality of the two 
expressions in Eq. (3) results from the definition of maturity 
in Eq. (1). 

In terms of practice, Eq. (3) implies that either the an¬ 
nual ring width or the bole radius must be modeled as a 
function of maturity to determine the dimensionless growth 
rate. It is somewhat easier to model the bole radius; and so 
we intend to determine the dimensionless growth rate using 


Rapidly grown trees have a lower heartwood content only 
when they are compared with more slowly grown trees of 
the same size. 13,14 When the age is the same, trees of large 
diameter appear to have a higher heartwood content. 13,14 Up 
to now, chronological age being considered the primary 
measure of the maturity of pine trees, the authors have had 
significant difficulty interpreting such observations. They 
appear to be contrary to findings according to which the 
proportion of heartwood increases with age, 6,8 " 11 and is 
negatively correlated with the growth rate of the tree. 6 ' 8,11,12 

Retarded maturation of trees grown on poor soil and in 
a cold climate 21-24 indicates that at a specified maturity, trees 
from regions and sites of inferior growth are older than 
trees that experienced more favorable growth conditions. In 
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Fig. 6. Normalized residual variation of the heartwood proportion as a 
function of the sapwood annual ring width (top) and as a function of 
cambium age (bottom) 


other words, at a specified chronological age trees from 
areas with favorable growth conditions are more mature 
than trees from areas of poor growth conditions. This ap¬ 
pears to explain the observation that large trees have a 
larger heartwood content than small trees at the same chro¬ 
nological age. 13,14 Werberg’s 13 data were gained with signifi¬ 
cant soil fertility variation, and Tamminen’s 14 data were 
obtained with significant climatic variation; both of these 
factors obviously induce a change in the time scale of the 
tree maturation process. 21 ~ 24 

The present approach to tree maturity appears to have at 
least one complication. The increment rate of size fluctuates 
owing to variations in environmental conditions, and so the 
relative increment rate of size does not decline steadily. 
This problem can be avoided by determining the increment 
rate of size over such an extended period of time that the 
fluctuations in the growth rate due to the variation in cli¬ 
mate are averaged. 

A suppressed tree may have a low growth rate over an 
extended period of time, and its growth can be significantly 



Fig. 7. Bole radius as a function of maturity. Eq. (4) is fitted to Dataset 
1 (squares). Logs of Dataset 2 ( circles) are somewhat larger at any 
specified maturity than logs of Dataset 1 



Fig. 8. Heartwood content as a function of dimensionless growth rate 
in Dataset 1 (diamonds) and Dataset 2 (triangles) 


enhanced if the suppressing trees are removed. Even if 
the rate of size increment were averaged to exclude fluctua¬ 
tions in climate, the inverse of the relative incremental rate 
of growth within such a suppressed tree may decrease with 
time. Such a behavior obviously is physical and so it must be 
accepted that maturation may be partly reversible. The sup¬ 
pressed tree actually may have been not far from death in its 
suppressed state, but it may have a significant amount of life 
ahead after removing the suppressing trees. 
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